# -*- coding: utf-8 -*-
# @Author  : longbhu
# @Time    : 2025/2/24 13:06
# @Function: 将一个完整的生态图，按照ecosystems_type.json划分成9大生态系统加2个水面。
# todo           灌丛归到森林生态系统，其他陆域不计算。
# todo           没有生态系统值的，不输出tif。生成tif后给出一个报告。
from calc_v1.utils import process_ecosystem_tif, extract_ecosystem_arrays

# 示例调用
tif_files = process_ecosystem_tif(
    r"F:\code\dev\calc-gep-regulate-cqc\devCalc\data\A003_生态系统分类栅格图_商洛-GCS_WGS_1984.tif"
)
print("\n生成的TIF文件列表：")
print(tif_files)

# 示例调用 extract_ecosystem_arrays
ecosystem_arrays = extract_ecosystem_arrays(
    r"F:\code\dev\calc-gep-regulate-cqc\devCalc\data\A003_生态系统分类栅格图_商洛-GCS_WGS_1984.tif"
)
print("\n各生态系统类型的numpy数组列表：")
for name, array in ecosystem_arrays.items():
    print(f"生态系统: {name}, 数组形状: {array.shape}")
